Particle-Hole Asymmetry in Doped Mott Insulators: Implications for Tunneling and 

Photoemission Spectroscopies 
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In a system with strong local repulsive interactions it should be more difficult to add an electron 
than to extract one. We make this idea precise by deriving various exact sum rules for the one- 
particle spectral function independent of the details of the Hamiltonian describing the system and of 
the nature of the ground state. We extend these results using a variational ansatz for the supercon- 
ducting ground state and low lying excitations. Our results shed light on the striking asymmetry in 
the tunneling spectra of high Tc superconductors and should also be useful in estimating the local 
doping variations in inhomogeneous materials. . 



With the discovery of high temperature superconduc- 
tivity in the cuprates there has been enormous inter- 
est in the properties of doped Mott insulators. One- 
particle spectroscopies like angle resolved photoemis- 
sion (ARPES) and scanning tunneling microscopy 
(STM) (3|,|Ui have played a major role in our understand- 
ing of these strongly correlated materials. In this pa- 
per we examine in detail sum rule constraints on the 
single-electron spectral function, focusing in particular 
on the striking asymmetry between occupied and unoccu- 
pied spectral weights in lightly doped Mott insulators and 
its doping dependence. Anderson Q has recently empha- 
sized that such an asymmetry is very unusual in conven- 
tional metallic systems, especially on low energy scales, 
and has suggested that it may be a characteristic sig- 
nature of "projection" into a low-energy subspace where 
strong local Coulomb repulsion makes double occupancy 
at a site energetically prohibitive. 

In this letter, we first make this idea precise by deriving 
several exact sum rules for the T = spectral function for 
both the occupied and unoccupied spectral weights in the 
low-energy subspace. These results are very general and 
do not depend either upon the details of the Hamiltonian 
or on any assumptions about the nature of the ground 
state or low-lying excitations. We show that these gen- 
eral results can be useful in several ways two ways. First, 
they help in quantifying the particle-hole asymmetry in 
tunneling spectroscopy: there is much more weight on 
the negative bias (occupied) side than on the positive 
bias (unoccupied) side [g- Second, our results can be 
used in present day STM experiments to estimate the 
local doping variations in inhomogeneous systems 2,]- 

We next separate out the coherent quasiparticle (QP) 
and the incoherent parts of the spectral function and pre- 
dict their variation with hole doping in the d-wave super- 
conducting state obtained upon doping a Mott insulator 
Hi This second set of results are variational in 

nature and require us to make assumptions about the 



ground state and low-lying QP excitations, and also to 
make the Gutzwiller approximation [lol llll | to obtain 
analytical results. These results are testable in photoe- 
mission and inverse-photoemission experiments. 

Exact Sum Rules 

Consider a system of electrons described by the Hamil- 
tonian H = K + U'^^ni^-Uji where K is the Kinetic 
energy operator which can be an arbitrary tight bind- 
ing Hamiltonian with terms of order t (nearest neighbor 
hopping). We can further add to _ff a random one-body 
potential which could make the system inhomogeneous, 
as well as other potential energy terms such as longer 
range Coulomb interactions. We will always work in the 
limit where U is much larger than all other energy scales. 

The one-electron spectral function is defined by 
A{r,r';uj) = ^ ImG'(r,r';aj + iO+) where G is the 
Green's function. We work in real space for two reasons: 
first, the no-double occupancy constraint is best writ- 
ten in this basis, and second, this allows us to describe 
spatially inhomogeneous systems which is important to 
discuss STM experiments on the cuprates. We use the 
T = spectral representation 



A{r,r';u;) = Y,[iO\clM)Hcra\0)S{uj + E^~Eo) 

m 

+ J2 [{0\Cra\m){m\cljO)SiLO - Era + £^0)] (1) 
m 

where |to)'s are exact many-body eigenstates with energy 
Efn with TO = the ground state, and uj is measured with 
respect to the chemical potential. 

The large U suppresses double-occupancy at each 
site, and its effects on the ground state and low-lying 
excitations are best described using the the projec- 
tion operator V = Yir (1 ~ Ur^Uri). We then make 
the well-known unitary transformation 12] exp(— zS'), 
such that exjp(iS)H exp{—iS) has no matrix ele- 
ments connecting states which differ in their double- 
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occupancy, to any given order in t/U . To leading or- 
der iS = -(l/t/)X]r,r',o-*rr' ('^rffCjcrCr'cr/lr'F-h-C.) + 

0{t/UY, where hr„ — l — n^.a- It is useful to incorporate 
this unitary transformation on the states, which is equiv- 
alent to transforming all operators. It then follows that 
all the low-energy states^ i.e., those in the so-called "lower 
Hubbard band" (LHB), are of the form exp(-iS')7'|$,„) 
where |<I>m)'s are wnprojected states. This characteriza- 
tion of LHB states will be crucial below. 

We now derive various exact sum rules without making 
any assumptions about the nature of the ground state or 
low-lying excitations. Some of these are very well known 
and shown only for completeness because we need to ref- 
erence them later. From eq. Q it is trivial to see the 
total spectral weight JJ*""" c?a;A(r, r'; — 1, while the 
occupied spectral weig 



/ (LjA{y,y'\uj) 

J —oo 



(2) 



Without assuming translational invariance, the local 
density of states (LDOS) probed by STM experiments 
7V(r; lS) — 274(r, r; w), with the factor of two coming from 
spin, is given by 



dLLjN{r; ut) — n(r) = 1 — x{r). 



(3) 



Here n(r) is the local electron density and x{r) the local 
hole doping, which for a translationally invariant system 
would be r- independent. This result simply says that 
there are (1 — x) occupied sites (per unit volume) from 
which one can remove an electron. 

For a translationally invariant system, we can Fourier 
transform (|2J) from (r— r') to k and obtain the well known 
result dujACkjLu) — n(k ) which has proved useful in 
analyzing ARPES data Summing over all k's and 

both spins one obtains dujN{uj) — 1 — x, which is Q 
for a uniform system. 

Next wc turn to sum rule constraints on the unoccu- 
pied side, which is the positive bias side in tunneling or 
that probed by inverse photoemission. It is trivial to 
derive sum rules for energy integration from (chemi- 
cal potential) to oo by subtracting the occupied spectral 
weights (0) or © from the total spectral weight of unity. 
But a much more physically meaningful result is obtained 
by focusing only on the low- energy states in the "lower 
Hubbard band" (LHB) by integrating over < a; < fi^, 
where the upper cut-off satisfies t <C JIl ^ U. This 
is implemented by restricting the sum over intermediate 
states in to LHB states |m) — exp(— z5)P|<i>m), as 
discussed above. 

We thus write the integrated low energy spectral 
weight on the unoccupied side as /g dujA{r,r';uj) = 

Em{'^o\'P5raV\^,n){^,n\V5l^V\^o)- Here we have 
found it convenient to move the unitary transforma- 
tion back onto the operators and have introduced 



the notation: c, o- = exp{iS)Crcr exp(— iS') and cj^ = 
exp(iS')4^exp(-iS'). We now use Eml*'")(^™l = 1' 
since the |$m)'s are wraprojected states, to obtain 



dujA{r, r'; uj) = ($o|^5r.^4'.^l*o). (4) 



must calculate 
VcLP = 



h -ct V 



^^'^ra'nraClcrV and 



is the 



To simplify this, we 
To order t/U we find: 

hermitian conjugate, where h^a = 1 — rir^a and a = —a. 

We thus obtain the sum rule for low energy spectral 
weight on the unoccupied side: 



dujN{r;uj)^2x{Y) + 2\{K{v))\/U (5) 



where {K{v)) = ($o|PER,.^Rr (cLcr. + h.c.) P|a>o). 
The first term in simply says that one can inject an 
electron into any of the x empty sites, with the factor 
of two for spin degeneracy. The second term gives an 
order (xt/U) correction since the injected electron can 
create a temporary double occupancy and then hop off 
to a neighboring empty site. We note that, in contrast to 
this, the corresponding result Q to extract an electron 
is exact to all orders in t/U. 

For a translationally invariant system, we may rewrite 
the above result as 



E 



dLuA(k,uj) = x+\{K)\/U. 



(6) 



There is another simple result that can be obtained in 
the translationally invariant case. First we simplify the 
right hand side of Q using the lowest order expressions 
for Vc)V and VcP. A straightforward calculation then 

shows that f^"- dujA{v,v';Lo) = (1 + .x)/2 - (cl,^Cra) + 
0(t/U). Fourier transforming to k-space and using stan- 
dard expressions of n(k), we find that the the total low- 
energy spectral weight is 



dcjyl(k, a;) = 



0{t/U) 



(7) 



for each k. Note that the deficit from unity comes from 
spectral weight in the "upper Hubbard band" which lies 
above Q^. 

Variational and Gutzwiller approximation results 



We emphasize that no approximations were made to 
obtain the above results, and we also made no assump- 
tions about the nature of the ground state or low-lying 
excitations. 

We now turn to translationally invariant systems and 
our goal is to obtain more detailed information about 
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the spectral function: to decompose it into its coher- 
ent and incoherent pieces and determine their doping 
dependence. We take the (variational) ground state to 
be a projected d-wave BCS state |0) = exp(zS')P|dBCS) 
which has given much insight into the phenomenology 
of the superconducting state of the high Tc cuprates 
P, "9]. Further the (variational) quasiparticle (QP) ex- 
citations llj above this ground state are described by 
I fed) — exp(z ^)^7LldBCS), where 7^ is the standard 
Bogoliubov QP operator. The QP's lead to the coherent 
part of A(k, uj), i.e., delta- functions in uj at T=0. Fin ally 
we make the Gutz wilier approximation (GA) 0, 0| 
which greatly simplifies calculations of matrix elements in 
this strongly interacting system and gives answers which 
are in good agreement with exact Monte Carlo results. 

Making a Gutzwiller approximation (GA) for the QP 
matrix elements in we obtain [l^ 

A{k,uj) = Z{k)ul6{Lj - Ek) 

+ Z(k)vlSiLu + Ek) + A„e(k, u) (8) 

where ttk, and i?k is standard BCS notation. The first 
two terms in (jSJare the coherent QP pieces with the same 
structure as in BCS theory, except that their spectral 
weight is suppressed by 

k' 

where ek is the dispersion corresponding to the bare ki- 
netic energy K in the Hamiltonian. We note that, as em- 
phasized in 1^, Z vanishes as one goes to the insulating 
state at a; = 0, and, in fact, the GA result © is in excel- 
lent quantitative agreement with the variational Monte 
Carlo results of Paramekanti et al. The sum over all 
states other than single QP's in Q leads to the incoher- 
ent part of the spectral function denoted by Ainc{k,uj). 
Although we cannot calculate its explicit form with the 
minimal set of assumptions we have made, its existence 
is necessarily demanded by exact sum rules, as shown 
below, which also put constraints on Ai„c- 

Our use of the GA to calculate matrix elements goes 
beyond previous applications of this approach, which 
have been by-and-large restricted to ground state ex- 
pectation values (see, however, the work of Laughlin 
who uses an approximation scheme closely related 
to GA for spectral functions). A non-trivial consistency 
check is provided by n(k) calculated within GA, which 
involves only an equal-time ground state correlation and 
does not depend on any assumptions about QP excited 
states. We find 7i(k) = Z(k)u^ + nsniooth(k) where 
nsmooth(k) = (1 - x)^/2{l + x) + 0{t/U) is a smooth 
function of k in the entire Brillouin zone. We omit the 
details of the (t/U) corrections here since they involve 



rather long expressions [ig. We note that n(k) implies 
that there is a jump discontinuity along the zone diago- 
nal whose magnitude is given precisely by JSJ including 
the (t/U) corrections. 

We now turn to sum rule constraints on Amc restrict- 
ing ourselves, for the most part, to leading order results 
in t/U; the next order corrections will be presented else- 
where 0. We begin by integrating ((HJ from —00 to 
and comparing with the GA result for n(k). We thus find 
that for each k 

£^ dc.A_(k,c.) ^ 1^:^ + 0{t/U). (10) 

Thus for each k there is non-zero incoherent spectral 
weight for w < 0, whose strength relative to the coher- 
ent weight Z(k)w^ on the occupied side, increases with 
underdoping (decreasing x). 

To find the incoherent spectral weight on the unoc- 
cupied side we substitute the GA spectral function ^ 
in the total low-energy spectral weight sum rule l|7|l. 
leading to J^^dcuA.nciKoj) = (1 - x)V[2(l + x)] + 
0{t/U) for each k. This together with ((Tn|l implies that 

P^'"- duoA,nc{K^) = 0{t/U). Given the non-negativity 
of spectral weight, we find that 

A„,(k,c^>0) = 0(^(7) <1. (11) 

The vanishing of Ainc for w > to zeroth order in {t/U) is 
at first sight quite surprising. Although there is very little 
spectral weight on the unoccupied side, as seen from ((HJ, 
whatever there is, is entirely dominated by the coherent 
piece when U ^ t. 

To gain more insight into this striking result we de- 
rive it in a completely different fashion, which also shows 
that it is not an artifact of the Gutzwiller approxima- 
tion. To zeroth order in (t/U) we can set exp(zS') = 1, 
and using the identity Vcl^P = T'cJ^, we find that 

VcljO) = Pc[^|dBCS) = constant x P-fljdBCS) = 
constant X |k, a). We thus find that the projected creation 
operator acting on the ground state gives precisely the co- 
herent QP state. Thus there is no incoherent weight in 
the electron creation (i.e., unoccupied) sector at least to 
zeroth order in {t/U). 

Finally, we can determine the explicit form of the order 
t/U incoherent unoccupied spectral weight, which follows 
from use of eqs. 10 and (TJ together with the GA result 
[3 for n(k). We find 

. ^ 2x(l - a:) o 

da,A_(k, ^) = ^7^^ E (12) 

We have also obtained sum rule constraints on the 
density of states contributions coming separately from 
the coherent and incoherent parts of the spectral func- 
tion, however we omit these here. In any case, in a k- 
integrated probe like tunneling it would be very hard to 
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separate out the coherent contribution from the incoher- 
ent one in experiments, unHke in ARPES where it seems 
possible to do so. 

Implications for Experiments 

We first discuss the impHcations of our results for tun- 
neling spectroscopy. The tunneling conductance in STM 
experiments is proportional to the local density of states 
G{r;eV) = constant x N{r;uj = eV) 18j where the con- 
stant of proportionality involves tunneling matrix ele- 
ments. Thus our result © shows that the (energy in- 
tegrated) positive bias conductance is small, of order x, 
while implies that the (integrated) negative bias con- 
ductance is large, of order unity. This provides a quali- 
tative explanation for the large asymmetry seen in STM 
experiments which show a superconducting gap structure 
superimposed on a sloping "background" which decreases 
going from negative to positive bias. Our results pre- 
dict how this asymmetry should grow with underdoping 
(decreasing x). This asymmetry is most strikingly seen 
in the highly underdoped ?7.on-superconducting cuprates 
such as NaxCai_xCu02Cl2 studied by Hanaguri et al. 
0. The nature of the "zero temperature pseudogap 
state" in such materials is an unsolved problem, and in 
this context it is very important to re-emphasize that 
our results Q and (j^l make no assumptions about the 
broken symmetry in the ground state or the nature of 
low- lying excitations. 

In order to get quantitative information from STM ex- 
periments we look at ratios in which the unknown tun- 
neling matrix elements cancel out. Taking the ratio of 
the total unoccupied low-energy spectral weig ht 10 to 
the total occupied spectral weight © we obtain 

j^^ <h:G{v-u:) 2x(r) 2\{K{v))\ 
f_^dLuGiv;Lu) [l-x{r)]^ U[l-x{r)Y ^ > 

The left hand side can now be estimated from STM data, 
provided one can make a reasonable choice of the positive 
and negative high energy cutoffs [l^, and then used to 
infer the local hole doping x{r) from the first term on the 
right hand side of H13|l . The second term of order {xt/U) 
gives an estimate of the approximately 10 % error made 
in estimating x. 

In the second part of the paper we derived results for 
the doping dependence of the coherent and incoherent 
parts of the spectral function. The predicted x depen- 
dence of the coherent weight Z of eq. © has already been 
observed in ARPES studies of nodal QP's (see ref. |^). 
The important prediction for the very small incoherent 
spectral weight on the unoccupied side pill2|) should 
be testable in future inverse photoemission experiments 
once their energy resolution is improved. 
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